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We consider a spin Hamiltonian describing d-d exchange interactions between localized spins d 
of a finite antiferromagnet as well as p-d interactions between a conducting hole (p) and localized 
spins. The spin Hamiltonian is solved numerically with use of Lanczos method of diagonalization. 
We conclude that p-d exchange interaction leads to localization of magnetic polarons. Quantum 
fluctuations of the antiferromagnet strengthen this effect and make the formation of polarons local- 
ized in one site possible even for weak p-d coupling. Total energy calculations, including the kinetic 
energy, do not change essentially the phase diagram of magnetic polarons formation. For parame- 
ters reasonable for high-T c superconductors either a polaron localized on one lattice cell or a small 
ferron can form. For reasonable values of the dielectric function and p-d coupling, the contributions 
of magnetic and phonon terms in the formation of a polaron in weakly doped high-T c materials are 
comparable. 
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I. INTRODUCTION 

The role of electron phonon and spin exchange in 
the formation of polarons, bi-polarons and stripes in 
CuC>2 based high-T c materials has been a matter of ex- 
tensive discussion from the very moment high-T c su- 
perconductors were discovered. Experimental observa- 
tions of the anomalous isotope effectErcl could indicate 
the role of phonon interactions in the formation of the 
superconducting state. n On the other hand, .-angularly 
resolved photoemissiorfl (ARPES-), transportcl and tun- 
neling microscopy measurementsu show d-symmetry of 
the order parameter indicating an important role of ex- 
change. Although, some authors underline the role of 
both mechanismsQ, the solution of full Hamiltonian is 
complicated and these contributions are usually calcu- 
lated separately. 

The concept of the ahonon polaron was introduced 
by PekarQ and FrohlichEna decades ago and was adopted 
to high-T c superconductors by Alexandrov and MottE2l. 
The energy of phonon coupling for high-T c is estimated 
to be of the order of a fraction of eV. In another 
work, AlexandrovEil proposed that the narrowing of d 
-electron band by the polaron effect and the formation 
of phonon bipolarons increases the critical temperature 
to T c w 100 K. However, because of the dominant con- 
tribution of (i-symmetry to the order parameter of high- 
Tc superconductors, exchange interactions are the most 
often considered ones recently. The modeling of such in- 
teractions is difficult because it is necessary to include 
the role of conducting (p) holes as well as localized (d) 
electrons. 

In the case of weak p-d couplings a linear response 
of the antiferromagnet (AF) can be assumed. Hence, 
a phenomenological model of magnetic susceptibility 



introduced by Millis, Monien, PinesB (MMP) can be 
used. 

Zhang and RiceEl, in turn, describe the case of strong 
hybridization between the 3c? states of Cu and 2p states 
of O. When the hybridization is much greater than the 
kinetic energy of the carriers, it is possible to limit the 
description to the ground singlet state, and in such a 
case a one band model can be used. This model (called 
t-J) is at present the most often used to describe the 
cuprate superconductors (for a review seeEil). Unfor- 
tunately, different techniques of approximations lead 
sometimes to contradictory solutipjas of t-J , as it is in 
the case of the stripe formatiorJl3'Ea. Also, the dynam- 
ics of spin polaron in the framework of t-J mqdek-was 
extensively investigated by many authorsOEJOEj' 2 3c3. 

However, the spin-fermion modc]il0S@0SS@, 
which treats separately the spin of p-carrier and d- 
localized spins, seems to be the most appropriate for 
the case of intermediate p-d coupling typical for high-T c 
superconductors. Most papers treat the p-d interaction 
as a Kondo exchange. Such an approach allows to sijow 
the tendency af-ftw>rates to stripe phase fermatiorcij or 
d-wave pairingEaEll Recently, Bala et al.c3£j extended 
the spin-fermion model starting from the three band 
Hubbard model. They introduced p-d coupling, which 
is much more complicated than Kondo one, i.e. some 
exchange couplings between Cu ion and two neighbor- 
ing oxygen atoms appear. They have considered-|tlie 
formation of magnetic polaron in such a modeEj'El 
They identified all .-the spectrum of energy .states of 
Zhang-Rice polaronEl in ARPUS experiments^!. They 
found a solution using self-consistent Born approxima- 
tion (SCBA). However, a comparison of SCBA and ex- 
act diagonalization techniques for spin-fermion model 
has not been carried out yet. 



The present paper reports for the hrst time on calcu- 
lations of a simple spin-fcrmion model with exact diag- 
onalization technique for spin Hamiltonian. It broadens 
the possible range of solutions. The spin-fermion model 
allows to investigate the cases of strong and intermedi- 
ate p-d coupling. Thus a comparison of the phonon and 
magnetic contributions to the formation of a polaron in 
weakly doped high-T c superconductors for a wide range 
of parameters is possible. The magnetic exchange in- 
teractions are calculated strictly while the phonon con- 
tributkjn_ia evaluated in the framework of the Frohlich 
theoryffiM 

In Section 2, we describe the Hamiltonian and the 
method of finding the ground and excited states of the 
system. In Section 3, the exchange spin interactions 
are analyzed. In Section 4, we construct and analyze 
phase diagram for the formation of different magnetic 
polarons. In Section 5, we compare the energies of 
phonon and magnetic polarons. A summary is given 
in Section 6. 



II. MODEL 

We study a system which consists of spins localized 
on the d-shells and mobile p-like holes. The localized 
spins are treated within the quantum antifcrromagnet 
model, i.e. none of the Neel sublattices are chosen a pri- 
ori. We consider finite one- (ID) and two- dimensional 
(2D) clusters as well as the infinite system. The latter 
is simulated using periodic boundary conditions (PBC) 
and an extrapolation to an infinite size. The spin in- 
teractions in the system and the kinetic energy of the p 
hole is described by the Hamiltonian: 

<i,j> i,ct,l3 

where Si is the effective spin of the nine core ci-electrons 
(3d 9 ), the summation was taken under the pairs of 
neighboring spins < i,j >, <r a p are the Pauli matri- 
ces of the carrier spin s, J p d is the Kondo parameter 
normalized to the volume of elementary cell, Jdd is the 
exchange integral between d-d spins, iy is the hopping 
matrix element to the neighbor, cf a , ci a are the oper- 
ators of creation and annihilation of a hole on the i-th 
site, with the projection of spin 1/2 on the quantization 
axis equal a respectively. We look for the eigenstates 
of the Hamiltonian (1) in the form of the product of 
the orbital wave vector, and the spin state vector of 
local and carrier spins X$: 



V-X s . (2) 

In the space of the single hole occupation at antiferro- 
magnet (AF) sites i = 1, N 

V = J2<P(i)ui, (3) 

i 

where Ui = |0, 0, 1, 0, 0...) is the occupation vector 
of the hole at i-th site and ip(i) is the corresponding 
probability amplitude. We use the variational method 
of solving the eigenproblem. We treat \P as the trial 
function, which is normalized by the condition: 

N 

X>«i a = i- ( 4 ) 

i 

For W given by Eq. (2) the diagonalization of the kinetic 
and exchange terms of Hamiltonian (1) can be done 
separately. 

H'= (*\H\V) =T + H S , (5) 

The kinetic energy, T, does not depend on the spin state 
Xg and can be easily evaluated from the probability 
amplitude of the carrier tp (i): 

r = *|ti|-2|ti| *>(*>0')> ( 6 ) 

<i,j> 

if only the hopping to the nearest neighbors (NN), t\ is 
considered, and z is the number of NN. The term z\t\\ 
has been added to have zero kinetic energy for fully 
delocalized states. Calculations of the exchange energy 
need diagonalization of the spin Hamiltonian: 

H s = 2 Jdd S * S J + X^WS* 8 - ( 7 ) 

<i,j> i 

Here s is the spin operator of the carrier spin and the 
p-d exchange integrals are given by 

j P d{i) = J P d \^{i)\ 2 ■ (8) 

We use Lanczos method of diagonalization of the spin 
Hamiltonian (7). Diagonalization of matrices for AF 
cluster with N < 20 is possible with the PC com- 
puter. For larger clusters the calculation time increases 
rapidly. 

A detailed study of quantum coupling between spins 
within the polaron in an AF medium is the main goal of 
the paper. We focus our analysis on the gain E ex of the 
spin exchange energy. E ex is the difference between the 
energy of the ground state of the unperturbed AF E^ d 



2 



(when J p d — 0) and the exchange energy of the whole 
system when the carrier is coupled to AF Edd + E p d, 
i.e. 



E, 



E dd 



(Edd + Epd). 



(9) 



In this way, we analyze the zero temperature case only. 
We found that for a weak p-d exchange the gain E ex is 
a half of the p-d coupling energy, E ex — —E p d/2. 

The Hamiltonian includes neither effects of external 
magnetic field nor any other terms which could lead to 
the symmetry breaking and to the formation of Neel 
sublattices. As a consequence, the mean values of the 
magnetic moments of each spin of this system are ex- 
actly equal zero. Because of this, we describe the spin 
structure of the antiferromagnetic polarons (AFP) and 
the mechanism of the magnetic polaron formation by 
means of correlators between spins of AF (S<Sj) and 
between the AF spin and that of an additional hole 



(sS,) = (XslsSilXs). 



(10) 



The assumption that the wave function can be de- 
scribed by Eq. 2 is simplified. In general, the total wave 
function is a linear combination of wave functions given 
by Eq. 2. In particular, because polarons localized on 
two Neel sublattices are equivalent, at least the combi- 
nation of the pair of such states should be considered. 
This problem is discussed in section IIIG. 

Hamiltonian (1) does not contain a phonon term. In 
general, the coupling of the carrier to the lattice polar- 
ization and the exchange interactions act together to 
form the polaron. The phonon term is treated in Sec- 
tion 5 within the simplest possible approach. It is cal- 
culated froro-Frohlich interaction with optical phonons. 
It is knownEj that in the range of parameters typical 
for high-T c superconductors small phonon polarons are 
formed. Their energy is given by: 



E, 



ph 



qD& 



(11) 



where go 



1/3 



is the Debye momentum, 



l . 



; e, £o are dynamic and static dielectric constants 
respectively, and e is the electron charge. 



III. SPIN STRUCTURE OF MAGNETIC 
POLARONS 



We examine various types of trial functions ip. De- 
pending on the size and shape of the carrier distribu- 
tion, \<p\ 2 (Fig. 1), as well as on the value of exchange 




icr 



o.oooo-_ 

o.o 0:1 6.2 
1/N 



-sV 



0.1 



FIG. 1: The exchange energy gain, E ex , obtained for a 
Gaussian trial function is presented as a function of polaron 
radius r' p , for ID open chain (crosses), for ID chain with 
(PBC) (open squares) and for a 4 x 4 square with PBC 
(open stars). Open triangles show E ex (r' p ) for a comb- like 
carrier trial function (see Eq. |lj) and ID AF chain of 16 
spins. In the inset the extrapolation of E ex to infinite AF 
ID cluster is shown by filled triangles, squares and circles 



for r' = 0.75ao, 



1.5ao and r' p = Sao respectively, where 



ao is a lattice constant. In the icons arrows show the corre- 
lators between the carrier spin and spins of AF for a small 
AFP (r' p = 0), medium AFP (r' p = 0.6a ) and large AFP 
(r' p — 1.4ao). For a large AFP, the arrow size is multiplied 
30 times for clarity. Solid curves show trial functions. 



constants (Fig. 2) various types of magnetic polarons 
with very different spin structures can be found. 

A. Large AFP. Its size r' p is much greater than 
the lattice constant ao and the antiferromagnetic cor- 
relation range £. Thus, it can be described within the 
theory of the linear response. For very large polarons 
the local magnetization._M(r) becomes proportional to 
the local carrier densityEj. 

B. Medium AFP. Its size r' p is comparable or 
smaller than £, but bigger than 0.5ao (see Fig. 1). The 
induced polarization is of staggered character, and the 
carrier interacts with few local spins of different orien- 
tation. 

C. Small AFP It is localized within a single elemen- 
tary cell. For very strong p-d coupling it corresponds to 
the Zhang-Rice polaron. 

D. Ferron. It forms when the p-d exchange is strong 
enough to break AF bonds and to induce a local mag- 
netic momentH. 

E. Comb like AFP It is the case when the carrier 
is distributed on one Neel sublattice only. The stag- 
gered polarization is the dominant response of the AF 
mediumSLJ. 

The clear difference between AFP of different sizes 
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FIG. 2: The exchange energy gain, E ex , from formation of 
different magnetic polarons as a function of J' pd . E ex for a 
small AFP were obtained for a ID chain with PBC (open 
circles) and for a 2D 4x4 square with PBC (filled circles). 
E ex for a trial function evenly distributed on all AF spins are 
shown by open and filled squares for a ID chain with PBC 
and 2D 4x4 with PBC respectively. Results for a trial 
function distributed on two neighbor AF spins and a ID 
chain with PBC are drawn by filled triangles. The dashed 
lines show quadratic and linear dependencies of E ex on J' pd . 
The dotted line emphasizes superquadratic dependence of 
Eex(J'pd)- The inset shows E ex (J pd ) for a small AFP. 



is seen in Fig. 1, where the energy gain, E ex , is plot- 
ted as a function of the polaron size r' p . Gaussian trial 



function \ f{i) 2 = Aexp 
malization factor and r, - 



(n-r ) 



where A is a nor- 



ro is the distance from central 
spin, was used to calculate almost all curves. The curve 
plotted by triangles, however, was calculated using an- 
other trial function, and it represents the formation of 
a comb-like AFP. It will be discussed in part III E. 

The calculations were performed for weak p-d cou- 
pling for ID and 2D AF clusters. The number of AF 
bonds has been assumed to be twice as big for 2D as for 
a ID cluster, where the number of ligands is two times 
smaller. To compare the results for ID and 2D clusters 
we introduce the following quantity, J' pd — J p d/ (2zJdd), 
which describes the ratio of p-d to d-d coupling nor- 
malized by the number of AF bonds, 2z. Ratio of the 
p-d to the d-d exchange energies defined this way is the 
same for ID and 2D AF clusters. Moreover, J' pd <Si 1 
describes the range of weak p-d coupling, while J' pd > 1 
that of strong p-d coupling, independent of the dimen- 
sion of the AF cluster. J' pd is set 1/4 in Fig. 1. Here 
the small AFP corresponds to r' p < 0.4a o , medium to 
0.4a o < r' p < a a , and large to r' p > a Q . The dependence 
of E ex on the chain length is shown in the inset. The 



linear dependence of the E ex on the inverse AF cluster 
size 1/N makes the extrapolation to an infinite chain 
possible. The extrapolated values of E ex are shown by 
full circles in Fig. 1. 

The results of calculations for small and medium AFP 
in the ID case are reasonable even without extrapola- 
tion to infinite clusters. Because of that, we assume that 
also for the 2D case the energy of small and medium 
polarons can be found with reasonable accuracy. Un- 
fortunately, because of computational limitations, we 
can not carry out full extrapolation for large AFP in 
2D clusters. 

A comparison of numerical results for open chains and 
chains with PBC shows that the use of PBC does not 
lead to an improvement of the convergence of numerical 
results with the cluster size. In part, this is an effect 
of the normalization condition given by (Eq. 4) and of 
specific properties of the dangling spins at the cluster 
border. 



A. Large AFP 

The numerical results, with the characteristic depen- 
dence of the E ex decreasing as l/r' p for ID AF, can 
be approximated by an analytical solution in which 
the magnetic susceptibility of AF is described by two 
phenomenological parameters. The analytical solution 
gives a better description of large AFP while the com- 
parison with the numerical results allows to explain the 
physical meaning of the phenomenological parameters. 

For simplicity in the analytical approximation, we use 
an exponential distribution of carrier density, |</?(r')| 2 = 
2^exp(— |r' — r |)/r p , where r is the center of the 

density distribution and r p = 2r 2 '. We assume that 
the susceptibility in wave vector space can be described 
by the phenomenological formula introduced in the 
MMP modelEl, x(q) = x qo i+(q-q )^ 2 ' where q is thc 
wave vector, qo — ir/ao and corresponds to the stag- 
gered magnetization. The phenomenological parame- 
ters Xqo an d £ describe the staggered magnetic sus- 
ceptibility and the spin correlation length, respectively. 
An effective field corresponding to (p(r') is iJ e //(r') = 



ig^Tr v ex P<H r ' - r |/r p ). Within the model of con- 
tinuous media, an effective field acting at a point r' 
results in magnetization M(r) = J dr'x(r — r')H e ff(r') 
at the point r. Such an approach takes into account the 
non-local effects of correlated systems but neglects the 
atomic structure of AF. For large AFP, i.e. for r p ^S> 
M(r) takes a simple form: 



M(r) = 



l6ir 2 gp, B £, 2 9o r 



exp 



r-ro 



(12) 
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Thus M(r) induced by a large AFP is smooth and de- 
scribed by the same spatial dependence as the carrier 
density. In the icon in Fig. 1, the correlators between 
the carrier spin, s and the consecutive local spins, Si, 
for a large AFP are presented. They also show that for 
a large AFP the induced magnetization is smooth and 
proportional to the carrier density. In that sense the 
analytical solution is equivalent to the numerical one. 
E ex — l/2E p d for a large AFP and it takes a simple 
form: 

E ex = (Jpdao)2 2 X ^ . (13) 

2567r 2 ( 5 /i B ) i 2 qlr p 

The absolute value of E ex and its dependence on 
the d-d exchange and on the correlation length cannot 
be analyzed within the discussed phenomenological ap- 
proach. However, an evaluation of E ex can be done by 
numerical calculations. If one attributes the correlation 
length to the cluster size, then the comparison of Eq. |l3| 
with the numerical result, showing that the energy gain 
is almost independent of the chain length (see inset to 
Fig. 1), leads to the conclusion that Xq (X £ 2 1 ' Jdd for a 
large ID polaron. We found out that the character of 
the obtained dependencies does not change if a Gaus- 
sian trial function was applied. 

In the 2D case, the results of numerical integration for 
large AFP in a 2D system are characterized by a faster 
decrease of E ex with the polaron size in comparison to 
the ID case. The energy scales with l/r p instead of 
l/ry 

To summarize, the large AFP is characterized by a 
smooth magnetization of the antiferromagnet, propor- 
tional to the carrier density. Because for a 2D po- 
laron E ex and the kinetic energies scale with the inverse 
square of the polaron size, two scenarios are possible: (i) 
no polaron can be formed when T > E ex , or, in the op- 
posite case, (ii) the minimum of the total energy occurs 
for the polaron radius tending to zero (see also Fig. 6). 
In that case, however, it is no longer a large polaron 
and a different approximation should be considered. 




FIG. 3: (a) The correlator, |(sSo)|, between the carrier spin, 
and AF spin on which the carrier is localized as a function of 
Jp d for ID chain with PBC (filled circles) and 2D 4x4 square 
with PBC (open stars). Icons show the correlators between 
the carrier spin, s, and consecutive d spins of AF chain for 
corresponding J' pd . (b) The z-component of carrier spin, s z , 
of central spin So z , as well as the sum of AF spins ^T Si Z 
are shown for an open chain of 16 spins as a function of J' pd . 
(c) The correlator between the spin on which the polaron is 
acting and the NN (filled circles) and the NNN spin (open 
squares) for a chain of 16 spins. The correlator between the 
spin on which the polaron is acting and the NN spin for a 
2D 4x4 square with PBC is presented by stars. 



merical calculations predict a smooth magnetization for 
r' p > a but a staggered one for r' p < ao- 

For a weak p-d coupling, E ex increases with the 
square of J' d (see filled triangles in Fig. 2). However, 
calculations show that for J' pd > 1 a more complex 
mechanism occurs. The AF bonds broke and a ferron- 
like polaron forms. The details are described in part 
HID. 



C. Small AFP 



B. Medium AFP 

As it is shown in Fig. 1, with a decrease of r' p , E ex 
does not follow the l/r' p dependence, but for r' p < a Q 
a sharp increase of E ex is seen. For ID the increase 
of E ex ~ 10, as compared to the extrapolated l/r' p de- 
pendence, and for 2D is even bigger. The rapid increase 
of E ex shows that the medium AFP has a tendency to 
localization. The decrease of r' p is also accompanied by 
a change of the character of induced magnetization in 
the AF medium. As it is shown by icons in Fig. 1, nu- 



When the polaron radius becomes comparable to the 
size of the elementary cell, E ex saturates. The plateau 
in Fig. 1 corresponds to the case when the polaron is 
localized in a single elementary cell, i.e. r' p < clq. The 
numerical results weakly depend on the size of AF clus- 
ter and on boundary conditions. The magnetization in- 
duced by the small AFP is similar to that of the medium 
AFP, as it is shown by icons in Fig. 1. 

A dependence of the numerically calculated E ex on 
the normalized ratio p-d to d-d exchange J' pd is shown 
by circles in Fig. 2. There is no considerable difference 
between ID and 2D cases. For a weak p-d coupling 
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E ex oc J' pd . The upper limit of this quadratic depen- 
dence corresponds to the case when the p-d coupling is 
equal to the d-d exchange. For higher J' pd the energy 
increases linearly. 

The correlator (sSo) between the carrier spin, s, and 
AF spin on which the carrier is localized, So, is shown 
in Fig. 3a. In the range of weak p-d couplings, the mag- 
netization of the local spin, as seen by the carrier spin, 
(sSo), increases linearly and saturates for the strong 
J' pd . This corresponds to quadratic and linear depen- 
dence of E ex versus J pd , respectively. For a strong J pd , 
the correlator (sSq) = —3/4 that means AF spin is 
compensated by the spin of a carrier. 

The icons in Fig. 3a show the correlators between the 
carrier spin s and the consecutive d spins of AF chain 
for different J' pd . It is seen that a small polaron induces 
a staggered magnetization within the AF correlation 
range, f. But the amplitude of (sSo) varies with J' pd . 
The increase of the correlator | (sSo) | is accompanied by 
a reduction of the correlators between the carrier spin 
and other d spins. Also, the correlator of the spin So 
with neighboring spins S;, (SoSi), decreases with J' pd 
increase. The dependence of correlator for NN and the 
next nearest neighbor (NNN) are shown in Fig. 3c. The 
correlators are not affected by the presence of a weak 
polaron, but for stronger J' pd they decrease as 1/ ' J' pd - 

The ground state of a system consisting of a carrier 
with s = 1/2 and an AF cluster of an even number of 
local spins is a spin doublet, for any J' pd . However, de- 
pending on the value of J pd , the magnetic moment can 
be transferred from the carrier to the system of local 
spins. For a non-interacting spin s, or for a very weak 
p-d coupling, the carrier spin has a magnetic moment 
while the total moment of localized spins is zero. In 
Fig. 3b the z-component of the carrier spin, s z , the local 
spin 5*o 2 and the sum of all local spins J^i $iz are plot- 
ted as a function of J' pd . To distinguish the z-direction, 
we put a small magnetic field which does not perturb 
the spin structure of polaron. With an increase of J' pd 
the magnetic moment is transferred from the carrier to 
the local spins. For a very strong p-d coupling, the mo- 
ment of local spins J2i &iz saturates at value 1/2. At 
the same time, the moment of carrier spin and of the 
local spin, vanish. The magnetic moments of the lo- 
cal spins, ^iSiz, are distributed in a staggered way 
among the whole cluster. In the limit of strong p-d cou- 
pling, J' pd 3> 1 the small polaron is equivalent to the 
Zhang- Rice (ZR) one. Here the spins, s and So, are 
compensated and form a local singlet state, not cou- 
pled to other local spins, S^ (see Fig. 3c and the icon in 
Fig. 3a.). On the other hand, the formation of a central 
pair of compensated spins is accompanied by a transfer 
of the magnetic moment to the neighboring local spins 
and formation of staggered M(r) around the central 



cell. Such an effect was described in details, for exam- 
ple ir£J. The reconstruction of AF bonds in the vicinity 
of the small AFP is driven by a gain of d-d exchange 
energy. It is of the order of J dd . It is a small contribu- 
tion to the whole energy of the ZR polaron but becomes 
important for J' pd = 1 . Moreover, the reconstruction of 
the AF environment should be taken into account when 
interactions between separate polarons are discussed. 

To summarize, the small AFP in the range of interme- 
diate p-d coupling (typical for high-T c superconductors) 
is not equivalent to the ZR polaron. The carrier spin is 
also partially correlated with AF spins in the vicinity 
of local singlet. 



D. Ferron 

As it is shown in Fig. 1, E ex for a weak p-d coupling 
decreases with an increase of the polaron size. For big 
J pd , however, some new effects occur. As it is shown 
in Fig. 2, for a weak polaron evenly localized on two 
d-spins antiferromagnetically correlated E ex is smaller 
by a factor ~ 80 than in the case of a small AFP. With 
an increase of J pd , instead of the expected saturation, 
a faster increase of E ex for such a polaron occurs. 

An even more pronounced effect is seen for a homo- 
geneously distributed carrier density. As it is shown by 
squares in Fig. 2, E ex for a weak p-d coupling is practi- 
cally zero but sharply increases for J' pd ~ 2. Moreover, 
a further increase of J' pd is accompanied with a step- 
like increase of E ex . Studies of the spin structure show 
that this critical behavior is caused by a breakdown of 
the AF bonds and the formation of a magnetic moment 
in the antiferromagnetic medium. Such a magnetic po- 
laroxuis. known as ferron, and it was described decades 
agct3E3. In the whole range of J pd values there is no 
considerable difference between the ID and 2D cases. 

The spin structure of a ferron localized on two local 
spins is complex. The correlator of the carrier spin s 
with the two local AF spins S p , (sS p ), is negative (for 
antiferromagnetic p-d coupling) for the whole range of 
J pd . For a weak J pd , (sS p ) increases linearly, then satu- 
rates for very large J' pd approaching —0.5. The sharper 
increase of E ex , which is seen for J' pd = 2.5 is caused 
by a spin- flip of the two local spins, S p . The system- 
atic increase of E ex , faster than the linear increase, 
for J' pd > 4 originates neither from the correlation 
(sS p ) nor from the correlation between two local spins 
(S p iS p2 ) but results from a systematic reconstruction 
of the AF environment of the ferron. The correlation of 
each of the S p with surrounding spins, Si, gradually de- 
creases. Simultaneously, there appears a reconstruction 
of d-d couplings among the surrounding spins. 
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FIG. 4: The correlator between the carrier spin and one of 
AF spin (sSi) (a), the total spin moment of the AFP po- 
laron, S 1 *"*, (b), and the correlator between two nearest AF 
spins (SiSi+i) (c) as a function of J' pd for a homogeneously 
distributed trial function. The case of a ID chain of 16 spins 
with PBC is presented by circles, the 2D 4x4 square with 
PBC is shown by open stars. 



More apparent data about the mechanism of the for- 
mation of a ferron can be found if one assumes a ho- 
mogeneous distribution of the carrier on the whole AF 
cluster with PBC. The data are shown in Fig. 4. The 
correlator between one of the AF spins and the carrier 
spin as a function of J' pd is presented in Fig. 4a. The 
step-like behavior of (sSj) is visible for both ID and 2D 
cases. For J' pd < 2.5 correlator (sSj) = and E ex —0. 
For J' pd > 2.5 E ex increases in accordance with the ob- 
served steps of (sSj). In Fig. 4b the total magnetic 
moment of the polaron, = S{ Z + s z is presented. 
It is also characterized by steps occurring for the same 
values of J' pd as those of (sSi) when J' pd is large. For 
J' pd < 9 and 2D case the total magnetic moment is 
equal to 1/2. For J' pd < 2.5 this moment is located at 
the carrier, for J' pd > 2.5 it is a result of the coupling of 
the electron spin (s 2 = 1/2) with Si Z — 1 AF spin. 
The correlator between two neighbor AF spins (SiS; +1 ) 
is presented in Fig. 4c. It changes from —0.4 for a pure 
AF to 0.25 (both spins parallel) with an increase in the 
strength of p-d coupling. At the same time, Sf" 1 ap- 
proaches value 7.5 (all bonds broken). It is worth to 
note that the step- like behavior of (SjSi+i) is not a 
consequence of breaking of the subsequent bonds, but 
it is associated with steps of the total magnetic moment 
of the system. 

To summarize, the ferron can form when homoge- 
neous, Curie- Weiss magnetization is induced. Our re- 
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FIG. 5: (a) The mean value of the correlator between the 
carrier spin s and the spins of the AF sublattice on which the 
polaron is acting Sa, \{sSa}\ (open circles), and the spins 
of other sublattice S_g, (sSb) (filled squares) as a function 
of J' pd for a comb-like trial function, (b) The z-component 
of the carrier spin, s z (filled squares), the central spin Sq z 
(open triangles), as well as the sum of AF spins ^\ Siz (filled 
circles) for an open chain of 16 spins as a function of J' pd . (c) 
The correlator between the spins from the same sublattice 
(SiSi+2) (filled circles) and from two neighbor sublattices 
|(SiSj+i)| (open squares) as a function of J pd . 



suits complete well the classical model of ferron. 



E. Comb-Like AFP 

A homogeneous distribution of the carrier on many 
AF spins makes the E ex on opposite spins cancel each 
other out. Thus we introduced a carrier trial function 
which is distributed on one of the Neel sublattices only 
as intuitively energetically favorable. The carrier dis- 
tribution is described as follows: 

| 2 = /(r,)cos 2 (I 4 ) 

where /(r^) is an envelope function, u suall y Gaussian. 
We call such a polaron comb-like AFPE£lE3. 

E ex versus J' pd for the comb-like polaron with a Gaus- 
sian envelope function is presented by triangles in Fig.l. 
For a small r' p , the comb-like AFP is equivalent to the 
small AFP. For a bigger r' , the gain E ex decreases, but 
it is still much greater than that for large AFP. E ex 
for a large comb-like AFP is only three times smaller 
than that of a small AFP. This stems from a decrease 
of quantum fluctuations (see part IIIF). 
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The comb-like AFP induces a staggered magnetiza- 
tion in the AF medium. This is shown in Fig. 5a, where 
the mean value of the correlator between the carrier spin 
s and the spins of the AF sublattice on which the po- 
laron is acting Sa, | (sS^) | , and the spins of other sub- 
lattice Sb, (sSb), are presented by circles and squares, 
respectively. The values of the correlators increase lin- 
early for small p-d coupling, and they saturate when 
J' pd > 1. The saturation value for (sS^) = 0.25 and it 
is independent of the number of spins. The saturation 
value of |(sSji)| > 0.25 and depends on the number of 
spins. The difference between the magnetization of the 
A and B sublattices gives the net magnetic moment of 
AF. 

A part of the magnetic moment which is transferred 
from the carrier to the AF medium depends on J pd , as 
it is shown in Fig. 5b. Here the change versus J' pd of z- 
components of the carrier spin, s z , of one of the spins 
on which the carrier acts, Sq z , and of the total spin, 
Sig, is presented. For small values of p-d coupling a 
quadratic increase of the moment transferred to AF is 
observed. Si Z saturates at a value smaller than 0.5. 
Simultaneously, the z-componcnt of the carrier spin de- 
creases from 0.5 for weak p-d coupling to the saturation 
value of 0.15 for an AF cluster consisting of 16 spins. 
However, s z + J2i Siz — 0-5 is independent on J pd . 

The presence of the comb-like AFP on one of the 
Neel sublattices breaks the translation symmetry and 
leads to the formation of two magnetized sublattices. 
In Fig. 5c the correlators between NN, |(S*Si+i)|, and 
NNN, (SiS i+2 ), are presented as a function of J' pd . 
While |(SiSj+i)| decreases for large J' pd , the correlator 
(SjS i+2 ) increases. Thus, a comb-like AFP increases 
the AF correlation radius. In the limit of very large J' pd 
and not too small sizes of the AF cluster, the comb- 
like AFP transforms the quantum antiferromagnet into 
a classical one. 

Although E ex for the comb-like AFP is large, the 
probability of its formation is small because it is char- 
acterized by a big kinetic energy. Moreover, the comb- 
like AFP has a tendency to localize into small AFP. 
This comes from the fact that the kinetic energy is al- 
most the same for both types of polarons (as long as 
only t\ is taken into account) while E ex due to quan- 
tum fluctuations increases three times with localization. 



F. Quantum effects 

In this section we summarize the results concern- 
ing the formation mechanisms of different magnetic po- 
larons. We pay particular attention to the appearance 
of quantum effects. In paragraphs A-E, five types of 
antiferromagnetic polarons were presented. The forma- 



tion of weak small, medium, and comb-like polarons is 
associated with the induction of staggered magnetiza- 
tion. In contrast, the formation of a strong ZR polaron 
is caused by the compensation of the AF spin by the 
carrier spin. The mechanism of ferron formation is the 
breaking of AF bonds, which is equivalent to the induc- 
tion of a homogeneous moment. 

Classically, the only possible mechanism of AF po- 
laron formation at zero temperature is the turning of 
the spins in the direction of the effective field induced 
by the carrier spir£3i23. However, this mechanism does 
not explain the formation of the remaining polarons. 
The induction of staggered magnetization is directly as- 
sociated with the quantum treatment of the AF. In a 
quantum approach the ground state of AF is a com- 
bination of Neel states, hence no sublattice is distin- 
guished. The simplest measure of AF quantum fluc- 
tuations is the difference between the ground and the 
first excited states, which are both composed of Neel 
states. In the presence of staggered magnetization the 
Neel sublattices are distinguishable and the damping of 
spin fluctuations occursEa. 

The mechanism of compensation responsible for the 
formation of the ZR polaron should be also considered 
in a quantum approach. Classically, the carrier cou- 
ples the AF spin and such interaction does not de- 
stroy the AF order. As we have shown in Fig. 3c, in 
the quantum treatment of the antiferromagnet the af- 
fected AF spin may cease to be correlated with neighbor 
spins. It means that only the quantum approach al- 
lows to model the destruction of antifecromagnetic or- 
der which is observed in experimentstiH. Despite the 
quantum character of the discussed above mechanisms, 
there are some particular situations where the polarons 
can be described in the classical limit, but the quan- 
tum corrections have to be taken into account. These 
corrections are the result of scaling of spin fluctuations 
with the number of spins and seem to be important 
for the antiferromagnetic sign of J' pd . By a system- 
atic analysis of the dependence of E ex on p-d cou- 
pling in the range of very large J' pdl we have found 
that E ex — ^2ij P d(i){sSi), where the summation is 
taken over the N spins Si within the polaron size (see 
Fig. 2). The correlator depends on N and is equal to 
(sSi) = (0.25 + 1/2 N). Thus E ex changes from 0.75 J' pd 
for a polaron localized on one spin to 0.25J pd for a 
polaron interacting with an infinite number of spins. 
Hence, it is possible to calculate the energy of strong 
small polarons as well as ferrons classically, but the ap- 
propriate scalar spin multiplication (0.25) has to be re- 
placed by the quantum factor. It approaches the clas- 
sical value when N — > oo, whereas for N — 1 it is three 
times bigger. The same quantum factor 3 is lost when 
the ZR polaron is delocalized on one of the Neel sub- 
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lattices forming the comb-like polaron. 

The importance of the sign of the p-d coupling is 
shown in the inset to Fig. 2. Here E ex for a small 
AFP with either a positive or negative sign of J' pd is 
presented. Different slopes of E ex for ferro- and antifcr- 
romagnetic coupling are observed. It is a consequence 
of the fact that the carrier and AF spins can not cou- 
ple ferromagnetically to more than 0.25 while for anti- 
ferromagnetic coupling (sSj) approaches —0.75. Thus, 
only for antiferromagnetic coupling the corrections as- 
sociated with the number of spins in the AF system 
are important. In the case of ferromagnetic coupling, 
the classical energy calculations for small polarons and 
fcrrons are correct. 



G. Polaron mass 

In our approach (see Eq.(2)) we assume a simple de- 
coupling of the carrier kinetic energy and the spin de- 
gree of freedom. For all the types of polarons discussed, 
the minimum energy corresponds to the localization of 
the maximum of the trial function, \l/(r), exactly at 
the position of the local spin ro. But another equiva- 
lent minima correspond to the localization at another 
equivalent spin positions. For PBC, which simulate the 
translation symmetry, all minima are exactly equiva- 
lent. As the consequence, the linear combination of the 
simplified solutions is expected to be a better solution of 
the problem. It shows that the carrier motion and the 
spin structure are really coupled. To study that prob- 
lem correctly one should, however, diagonalize N times 
bigger matrixes. Anyway, our simplified approach gives 
a reasonable picture of the spin structure of a magnetic 
polaron under assumption of its static nature, but it 
does not allow to study the real correlation of the car- 
rier motion and spin dynamics and does not bring any 
solid answer on the problem of the polaron localization. 
However, a hint about the effective polaron mass and 
the reduction of the kinetic energy can be found within 
our simple approach through the analysis of the change 
of the spin structure which accompanies the transfer 
of the polaron to the NN. When the decoupling of the 
motion and spin is assumed (Eq. (2)), the kinetic en- 
ergy is described by the parameter t\ , according to Eq. 
(6). Here, only the transfer without any spin flip is con- 
sidered. Our solutions show, however, that the trans- 
fer is accompanied with a complex change of the spin 
structure. The spin structure of the surrounding local 
spins is very different when the carrier is localized on 
the first, A, or on the neighboring site, B. A posteriori, 
we are able to claim that the matrix elements describ- 
ing the kinetic energy are reduced by a factor resulting 
from the product of the spin states on the site A and 



B (Xg\Xg), respectively. In other words, to consider 
the effects which are omitted due to our assumption 
(Eq.(2)), one has to reduce the value of the parameter 
of the kinetic energy. Instead of t\ as effective kinetic 
energy, t\ = t\(Xg\Xg) should be used. Our numeri- 
cal study shows that for weak p-d coupling, J p d <C Jdd, 
the factor {Xg\Xg) is close to unity. For J p d w Jdd the 
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value of (Xg\Xg) decreases, and for J p d > Jdd it satu- 
rates at the value in the range 0.1 - 0.2, depending on 
the type of polaron. It shows that the effective mass of 
polaron dressed by the spin polarization is considerably 
bigger as compared to free polaron. 



IV. PHASE DIAGRAM OF MAGNETIC 
POLARONS 

Until now we have considered the E ex from mag- 
netic polarons formation and we have described the 
spin structure of various polarons. However, the to- 
tal energy (exchange, phonon and kinetic one) should 
be taken into account to specify which type of polaron 
forms. In this section we consider the formation of a 
pure magnetic polaron. 
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FIG. 6: Total energy, E to t, as a function of polaron radius 
r' p for J' pd — 1.875 and effective masses of m e (triangles), 
2m e (circles), 25m e (squares). 



In Fig. 6 the total energy (sum of kinetic and ex- 
change energies) as a function of polaron radius calcu- 
lated for a few chosen effective masses is shown for a 
2D square cluster with PBC. It can be seen that for 
m e f f = m e no AFP can form. For a large m e f / (which 
scales as 1/ti), the only polaron which can exist is the 
small one. The formation of a ferron is not possible 
because J' pd is too small. The medium and large po- 
larons tend to localize in one elementary cell, indepen- 
dently of the values of m e ff and t\. The tendency to 
localization is associated with a rapid reduction of E ex 
with increasing r' p , which stems from quantum fluctu- 
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ations of the AF medium (see part IIIF). We found 
that for J' pd — 1.875 magnetic polarons can form only 
if m e ff > 1.6m e , what corresponds to t\ — 2 Jdd- For 
weaker p-d couplings the small AFP is also the only 
polaron which can form but its m e ff should be much 
greater. 

In Fig. 7 the phase diagram of the formation of differ- 
ent magnetic polarons is shown. The area of parameters 
ti and J' pd for which the small AFP forms is marked by 
a solid line. 




FIG. 7: Phase diagram of magnetic polaron formation. In 
the inset Etot as a function of ti for J' pd = 4.4 is shown. In 
the inset dashed line presents Etot{ti) for two-spin ferron, 
the solid line shows Etot(ti) for small AFP. Vertical lines 
separate the areas of different polarons formation 

For higher J pd , the small ferrons can also form. In the 
inset to Fig. 7 the total energies of a small AFP and a 
small ferron are compared for J' pd — 4.4. E ex for a small 
AFP is higher than for the two spin ferron, and it is 
constant. The kinetic energy of two-spin ferron (polaron 
localized on two neighboring spins) is smaller than that 
of a small AFP and scales linearly with t\. Thus the 
total energies of both polarons are linear functions of 
t\ but with different slopes. In particular the boundary 
between small AFP and two spin ferron for J' pd = 4.4 
crosses for t\ = bAJdd, what is equivalent to m e // > 
0.6m e . For ti, which is greater than 6.25Jdd any AFP is 
not formed. In general, the small AFP forms for m e f t > 
m e , while the two-spin ferron appears for m e tf < m e . 

Similarly, we found the range of parameters for which 
other types of polarons can also form. In Fig. 7 the 
ranges where small AFP, two, four-spin ferrons can oc- 
cur are shaded by solid, dotted and dashed lines, re- 
spectively. The appearance of a three-spin ferron is in- 
dicated by rectangles. Larger ferrons appear for much 
stronger p-d couplings, and we omit their description 
because they would exist for nonphysical parameters. 

The presented results were calculated for a 2D square 



with PBC. To estimate the error caused by small size 
of antiferromagnetic cluster, we investigated also such 
clusters as: a 4 x 4 square without PBC, a 4 x 4 par- 
allelogram, and 2 x N systems where 2 < N < 10. We 
found that the energy of polarons with small size de- 
pends on the nearest vicinity only. It is weakly depen- 
dent on both N and boundary conditions. For different 
2D clusters the formation range of subsequent polarons 
can be shifted about 10%-15%. 

The investigation of polaron formation in a wide 
range of parameters is exceptionally valuable when we 
want to classify what type of polaron could be formed in 
different materials. For example,.-^ doped EuTe where 
J pd 0.125 eV and Jdd ~ 2 meVl2a the p-d coupling on 
one AF bonds J' pd ~ 8. At the same time, ti ps l5Jdd- 
For such parameters our model predicts the ferrons 
mation (see Fig. 7) in accord with experiments (seeE2 
and references therein). 

In high-T c superconductors which are based on CuC-2 
layers, the d-d exchange coupling ranges from 0.055 eV 
to 0.085 eV, as determined from the Neel temperature 
of antiferromagnetic precursors of these materialsE3. 
Jpd ~ 1 — 3 &Vj was, determined from the band structure 
calculationsacjOcd. Thus, the ratio of p-d to d-d ex- 
change energies is between 1.5 and 4 for Jdd — 0.075 eV. 
m e ff is much more difficult to determine because the 
mass observed in experiments (e.g. in ARPES) is dec- 
orated by phonon and exchange interactions. Since-the 
decorated m e ff observed in ARPES experimentsCJ is 
about 4m e and m e // in our model is dressed only by 
phonons, the m e // which should be considered in the 
phase diagram is smaller than 4m e , which is equivalent 
to ti > 0.8J dd . 

We conclude on the basis of our phase diagram that 
for parameters typical for weakly doped high-T c super- 
conductors mainly the small AFP would form. There 
may be also some competition between the formation 
of small AFP and two-spin ferrons provided that the 
effective mass dressed in phonons is a bit lower than 
the electron mass. Large AFP have energies of few or- 
ders of magnitude lower than those of small AFP and 
ferrons. 



V. THE CONTRIBUTION OF PHONON AND 
MAGNETIC POLARON 

Now we will compare the phonon and magnetic con- 
tributions to the formation of polarons in weakly doped 
high-T c superconductors. We have shown that in the 
2D case E ex from the formation of large magnetic po- 
larons is close to zero and that they are unlikely to form 
(Fig. 6 and Fig. 7). Mott and AlexandrovM suggest that 
in high-T c superconductors only small phonon polarons 
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FIG. 8: E ex from the formation of a small AFP as a function 
of Jp d (a). The phonon energy, E p h, versus effective dielectric 
constant k (b). The absolute value of the total energy for a 
magnetic polaron (circles) and for a phonon polaron (solid 
line) as a function of J' pd (c). Here, m e ff = 2.2m e and k = 5. 
The vertical dashed line shows the value of J' pd where the 
magnetic and phonon polaron energies are equal. 



In order to obtain the total energy of magnetic and 
phonon polaron separately, E to t, we have evaluated 
the energies for a chosen set of parameters: m e // = 
2.2m e and k = 5. This mass, according to_ARPES 
experimental and band structure calculations^], is rea- 
sonable for high-Tc superconductors, k is taken from 
theory of phonon polarorcj. In Fig. 8c, we show the 
values of E to t for phonon as well as magnetic polarons 
versus J' pd . For the phonon polaron E to t is constant 
and equal to 0.15 eV since the phonon as well as kinetic 
energies are independent of J' pd . The phonon polaron 
dominates in the range of weak p-d couplings. In con- 
trast, the magnetic polaron (circles in Fig. 8c) dominates 
in the range of strong J' pd where its energy increases lin- 
early with J' pd . In the range of J' pd — 1.5 — 4, which is 
typical for weakly doped high-T c superconductors, the 
energies of both polarons are comparable or magnetic 
polaron energy slightly overcomes the energy of phonao 
one. This result agrees well with a very recent papeiEa 
which stresses the importance of phonon interactions in 
high-T c superconductors. 

We conclude that within our model the phonon and 
magnetic polaron energies are comparable and that 
both contributions should be included in the total po- 
laron energy in high-T c superconductors. 



VI. SUMMARY 



form (see Eq. (12)). Thus we will not consider large 
polarons here, neither of magnetic nor phonon origin. 

In Fig. 8a, E ex due to the formation of a small mag- 
netic polaron in a 2D AF square with PBC is shown as 
a function of J' pd . In Fig. 8b, the energy gain due to 
the formation of phonon polaron is presented as a func- 
tion of effective dielectric constant n. Typical values of 
k are between 5 and 20. The kinetic energy for both 
small polarons is the same, Ek = 4|ti|, and is omitted 
for clarity. 

Horizontal and vertical dotted lines drawn in Figs. 
8a and 8b indicate parameters J' pd — 1 and k — 10 for 
which energy gain E ex — E p h — 0.35 eV is the same 
for both phonon and magnetic polarons. The energy 
gain due to the formation of phonon polarons is greater 
than that for magnetic polarons for k < 10 and J' pd < 1, 
while k > 10 and J' pd > 1 the magnetic polaron dom- 
inates. One can see that for the range of parameters, 
which seem to be typical for high-T c superconductors, 
the energy gain due to the formation of phonon and 
magnetic polarons is comparable and of the order of 
fraction of eV. 



Using the spin-fcrmion model and treating all spins 
in quantum approach, we found that depending on 
the material parameters five various types of magnetic 
polarons in antiferromagnetic medium can be distin- 
guished. 

We showed that in the range of parameters typical 
for weakly doped high-T c superconductors the contribu- 
tions of phonon and magnetic interactions to the forma- 
tion of a polaron are comparable and hybrid magnetic- 
phonon polarons form. 

In weakly doped high-T c superconductors small mag- 
netic polarons as well as small ferrons can form. Our nu- 
merical study allows to find a continuous evolution from 
Zhang-Rice approach, when p-d coupling is assumed to 
be much stronger than d-d coupling, via important for 
high-r c materials intermediate range, to linear response 
approach suitable for weak p-d coupling. 

The here introduced comb-like polarons cannot form 
when the band is empty (i.e. in undoped material) be- 
cause of the cost of high kinetic energy. However, they 
can play a fundamental role in dope d high-T c materials 
when the band is partially filledc3c3. 



11 



ACKNOWLEDGEMENTS 



Valuable discussion with H. Przybylinska. Work sup- 
ported by the KBN grant 2 P03B 007 16. 



1 M. K. Crawford, M. N. Kunchur, W. E. Farneth, E. M. 
McCaron, and S. J. Poon, Phys. Rev. B 41, 282 (1990). 

2 J. P. Franck, J. Jung, M. A.-K. Mohamed, S. Gygax, and 
G. I. Spoule, Phys. Rev. B 44, 5318 (1991). 

3 M. R. Norman, H. Ding, M. Randeria, J. C. Cam- 
puzano, T. Yokoya, T. Takahashi, T. Mochiku, K. Kad- 
owaki, P. Guptasarma, and D. G. Hinks, Nature 392, 157 
(2000). 

4 K. Karpihska, M. Z. Cieplak, S. Guha, A. Malinowski, 
T. Skoskiewicz, W. Plesiewicz, M. Berkowski, B. Boyce, 
T. R. Lemberger, and P. Lindenfeld, Phys. Rev. Lett. 84, 
155 (2000). 

5 B. G. Levi, Phys. Today 17, 53 (2000). 

6 C. Castellani, C. D. Castro, and M. Grilli, Phys. Rev. 
Lett. 75, 4650 (1995). 

7 S. I. Pekar, Zh. Eksp. Teor. Fiz. 16, 335 (1951). 

8 H. Frohlich, Phys. Rev. 79, 845 (1954). 

9 H. Frohlich, Adv. Phys. 3, 325 (1954). 

10 A. S. Alexandrov and N. Mott, Polarons and Bipolarons 
(World Scientific, Singapore, 1995). 

11 A. S. Alexandrov, Phys. Rev. B 38, 925 (1988). 

12 A. Millis, H. Monien, and D. Pines, Phys. Rev. B 43, 258 
(1991). 

13 F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 

(1988) . 

14 E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

15 S. R. White and D. J. Scalapino, Phys. Rev. Lett. 84, 
3021 (2000). 

16 C. Hellberg and E. Manousakis, Phys. Rev. Lett. 84, 3022 
(2000). 

17 S. Schmitt-Rink, C. Varma, and A. E. Ruckenstein, Phys. 
Rev. Lett. 60, 2793 (1988). 

18 C. L. Kane, P. Lee, and N. Read, Phys. Rev. B 39, 6880 

(1989) . 

19 V. Elser, D. A. Huse, B. I. Shraiman, and E. D. Siggia, 
Phys. Rev. B 41, 6715 (1990). 

20 A. Ramsak and P. Horsch, Phys. Rev. B 48, 10559 (1993). 

21 G. F. Reiter, Phys. Rev. B 49, 1536 (1994). 

22 Z. Liu and E. Manousakis, Phys. Rev. B 45, 2425 (1992). 

23 C. Buhler, S. Yunoki, and A. Moreo, Phys. Rev. Lett. 
84, 2690 (2000). 

24 B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Phys. Rev. 
B 52, 5563 (1995). 

25 E. M. Hankiewicz, R. Buczko, and Z. Wilamowski, Phys- 



ica B 284-288, 437 (2000). 

26 J. R. Schrieffer, J. Low Temp. Phys. 99, 397 (1995). 

27 A. Ramsak and P. Prelovsek, Phys. Rev. B 42, 10415 
(1990). 

28 M. Hamada and H. Shimahara, Phys. Rev. B 51, 3027 
(1995). 

29 J. Bala, A. M. Oles, and J. Zaanen, Phys. Rev. B 54, 
10161 (1996). 

30 J. Bala and A. M. Oles, Phys. Rev. B 61, 6907 (2000). 

31 J. J. M. Pothuizen, R. Eder, N. T. Hien, M. Matoba, 
A. A. Menovsky, and G. A. Sawatzky, Phys. Rev. B 78, 
717 (1997). 

32 T. Kasuya and A. Yanase, Rev. Mod. Phys. 40, 684 
(1968). 

33 E. M. Hankiewicz, R. Buczko, and Z. Wilamowski, Acta 
Phys. Polonica A 97, 185 (2000). 

34 M. Laukamp, G. B. Martins, C. Gazza, A. L. Malvezzi, 
E. Dagotto, P. M. Hansenand, and A. C. L. J. Riera, 
Phys. Rev. B 57, 10755 (1998). 

35 V. D. Lakhno and E. Nagaev, Fiz. Tverdogo Tela 18, 
3429 (1976). 

36 Z. Wilamowski and E. M. Hankiewicz, Acta Phys. 
Polonica A 97, 403 (2000). 

37 R. J. Birgeneau, D. R. Gabbe, H. P. Jenssen, M. A. Kast- 
ner, P. J. Picone, T. R. Thurston, G. Shiraneand, Y. En- 
doha, M. Sato, K. Yamada, et al., Phys. Rev. B 38, 6614 
(1988). 

38 P. Wachter, in Handbook on the Physics and Chemistry 
of the Rare Earths (1979), p. 507, eds. K. A. Geschneider, 
Jr. and L. Eyring. 

39 N. M. Plakida, High Temperature Superconductivity 
(Springer Verlag, Berlin, 1995). 

40 M. S. Hybertsen, E. B. Stechel, M. Schluter, and D. Jen- 
nison, Phys. Rev. B 41, 11068 (1990). 

41 H. Eskes, L. H. Tjeng, and G. Sawatzky, Phys. Rev. B 
41, 288 (1990). 

42 C. X. Chen, H. B. Schiittler, and A. J. Frcdro, Phys. Rev. 
B 41, 2581 (1990). 

43 G. Mante, R. Claessen, T. Buslaps, S. Harmand, 
R. Manzke, and M. S. J. Fink, Z. Phys. B 80, 181 (1990). 

44 A. Lanzara, P. V. Bogdanov, X. J. Zhou, S. A. Kellar, 
D. L. Feng, E. D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, 
K. Kishio, et al., Nature 412, 510 (2001). 



12 



